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The £> e — e~ elastic scattering cross-section was measured with a CsI(Tl) scintillating crystal 
array having a total mass of 187 kg. The detector was exposed to an average reactor v £ flux of 
6.4 x 10 12 cm _2 s _1 at the Kuo-Sheng Nuclear Power Station. The experimental design, concep- 
tual merits, detector hardware, data analysis and background understanding of the experiment are 
presented. Using 29882/7369 kg-days of Reactor ON/OFF data, the Standard Model (SM) elec- 
troweak interaction was probed at the squared 4-momentum transfer range of Q 2 ~ 3 x 10" 6 GeV 2 . 
The ratio of experimental to SM cross-sections of f = [1.08 ± 0.21(stat) ± 0.16(sj/s)] was mea- 
sured. Constraints on the electroweak parameters (<7v,<?a) were placed, corresponding to a weak 
mixing angle measurement of sin 2 6\v = 0.251 ± 0.03l(stat) ± 0.024(sys). Destructive interference 
in the SM v c — e process was verified. Bounds on anomalous neutrino electromagnetic properties 
were placed: neutrino magnetic moment at fj,p e < 2.2 x 10 _10 /xb and the neutrino charge radius at 
-2.1 x 10~ 32 cm 2 < (r| e ) < 3.3 x 10" 32 cm 2 , both at 90% confidence level. 

PACS numbers: 14.60.Lm, 13.15.+g, 25.30.Pt. 



I. INTRODUCTION 

The compelling evidence of neutrino oscillations from 
the solar, atmospheric as well as long baseline acceler- 
ator and reactor neutrino measurements implies finite 
neutrino masses and mixings [l]]. Their physical origin 
and experimental consequences are not fully understood. 
Experimental studies on the neutrino properties and in- 
teractions are crucial because they can shed light to these 
fundamental questions and may provide hints or con- 
straints to models on new physics. 

We report a study of neutrino-electron scattering us- 
ing reactor neutrinos at the Kuo-Sheng Nuclear Power 
Station with a CsI(Tl) scintillating crystal array. The 
cross-section formulae are summarized in Scction[n] The 
conceptual design, hardware construction and perfor- 
mance are presented in Section IIII1 followed by discus- 
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sions on event reconstruction, background understanding 
and suppression, as well as experimental systematic ef- 
fects. Section IVTII shows results on the Standard Model 
(SM) electroweak physics @ as well as constraints on 
possible neutrino electromagnetic interactions. 

II. NEUTRINO-ELECTRON SCATTERING 

Neutrino-electron scattering has been studied with sev- 
eral generations of experiments at the accelerator using 
mostly muon-neutrinos v^v^) 0, Bj|. It is a pure lep- 
tonic process and therefore provides a clean test to SM. 
The typical squared 4-momentum transfer was Q 2 ~ 
1CU 2 GeV 2 and the electroweak angle sin 2 #w was probed 
to an accuracy of ±3.6%. 

Using electron-neutrinos as probe, the interaction 

u e (v e ) + e" -> v c {v c ) + e" (1) 

has been studied at medium energy accelerators [!, 
as well as at the power reactors pHHj]. It is also an 
important channel in the detection of solar neutrinos 
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TABLE I: Summary of published v e — arid v c — e scattering cross-section and sin Ow measurements. Unavailable entries are 
denoted by "N/A". 



Experiment 




E„ (McV) 


T (MeV) 


Events [14] 


Published Cross-Section 


sin 2 #w 


Accelerator v c : 














LAMPF [5] 




7 < E„ < 50 


7-50 


236 


[10.0 ± 1.5 ±0.9] ■ E v 


0.249 ± 0.063 










xl0~ 45 cm 2 




LSND [6] 




20 < E„ < 50 


20-50 


191 


[10.1 ±1.1 ±1.0] - Ev 


0.248 ± 0.051 










xlQ- 45 cm 2 




Reactor i/ e : 














Savannah River 


















1.5 < E„ < 8.0 


1.5-3.0 


381 


[0.87 ±0.25] -a V -A 


f u.zyznu.uo 
} N/A 

o.22±8:| 


Original [7] 


{ 


3.0 < E v < 8.0 

1.0 <. -Cj^ <, O.U 


3.0-4.5 

1.0 — O.U 


77 

N/A 


[1.70 ± 0.44] • a V -A 

[1.O0 it U.4J • CTSM 


T Re-analysis [131 
Krasnoyarsk tgf* 


{ 


3.0 < E„ < 8.0 
3.2 < E v < 8.0 


3.0-4.5 
3.2-5.2 


N/A 
N/A 


[2.0 ± 0.5] ■ crsM 
[4.5 ± 2.4] 
x 10~ 46 cm 2 /fission 


Rovno [9] 




0.6 < E v < 8.0 


0.6-2.0 


41 


[1.26 ±0.62] 
xl0~ 44 cm 2 /fission 


N/A 














MUNU [10] 




0.7 < E„ < 8.0 


0.7-2.0 


68 


[1.07 ± 0.34] events/day 


N/A 


TEXONO (This 


Work) 


3.0 < E„ < 8.0 


3.0-8.0 


414±80±61 


[1.08 ±0.21 ±0.16] -crsM 


0.251 ±0.031 ±0.024 



Electroweak Parameters 





II- 



CC 




FIG. 1: Interactions of v c with electron via the SM-allowed 
charged current (CC) and neutral current (NC) channels. 
There is in addition interference effect between them. 



The SM differential cross-section in the laboratory 
frame for v^v^—e elastic scattering, where only NC is 
involved, is given by @,0|; 



— (We) 



SM 



—7. [ (9v ± 9 a) 

Z7T 

(gv T9a) 2 (i - -g- 



(g 2 v-g A ) 



El 



(2) 



where the SM v c — c scattering cross-section was used to 
extract neutrino oscillation parameters. This process is 
among the few of the SM interactions which proceed via 
charged current (CC), neutral current (NC) as well as 
their interference (Int) (HI, as illustrated schematically 
in Figure [T] The interference effect in v c — c scattering is 
the origin of matter oscillation of solar neutrinos in the 
interior of the Sun [l[ . 

The experimental results on v c — and v c — e scattering 
are summarized in Table HJ Neutrino-electron scattering 
was first observed with reactor reactors in the Savannah 
River experiment 0- Re-analysis of the data by a later 
work [131 ] with improved input on the reactor neutrino 
spectra and electroweak parameters gave cross-sections 
which were about 2<r higher than the SM values. The 
discrepancies were interpreted as hints of anomalous neu- 
trino interactions. Other subsequent experiments [8l-[l0| 
focused on the searches of neutrino magnetic moments 
at low recoil energy such that their sensitivities to SM 
physics were limited. 



where Gf is the Fermi coupling constant, T is the ki- 
netic energy of the recoil electron, E v is the incident 
neutrino energy and gv, gA are, respectively, the vector 
and axial-vector coupling constants. The upper(lower) 
sign refers to the interactions with v^T/p). For u e (D c )—c 
scattering, all CC, NC and Int are involved [l2|, and the 
cross-section can be obtained by making the replacement 
gv.A (gv.A ± !)• In the case of v c — c which is relevant 
for reactor neutrinos, 
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The SM assignments to the coupling constants are: 
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where sin 2 #w is the weak mixing angle. The SM differ- 
ential cross-section expressed in terms of sin 2 #vv is ac- 
cordingly: 
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The observables in an experiment are the event rates 
(Rexpt)- The SM predicted rate, expressed in unit of 
kg _1 day _1 , can be written as 



Rsm{v) = p e 



rd ° r ]sM 41 dE v dT , (6) 



dE„ 



where p e is the electron number density per kg of target 
mass, and d<f>u/ dE v denotes the neutrino spectrum. 

Results of this work are reported in several schemes 
using R expt ■ Firstly, the cross-section ratio 



Rsm{v) 



(7) 



can be used to probe new physics in a model- independent 
way. Alternatively, taking SM electroweak interactions 
but allowing the parameters to assume any values, the 
allowed ranges of (gv>9A) as well as sin 2 6>w can be de- 
rived from R eX pt, following Eqs.EMSl respectively. 

To study the interference effects, the measured rate 
can be expressed as 



R 



expt 



RCC + RNC + V ■ Rlr, 



(8) 



The CC— NC interference for v e (p e )— e is destructive in 
SM, or equivalently r){SM) = — 1. Possible deviations in 
the sign and magnitude of the interference effects (rj) can 
be probed. 

It follows from Eqs. [S$4H] and the analogous formu- 
lae for v B — e that, under realistic experimental configu- 
rations, the projected accuracies on sin 2 #w (denoted by 
A[sin 2 #w]) are related to the experimental uncertainties 
in £ (denoted by A[£]) by: 



A[sin 2 6» w ] ~ { 



0.15 
0.35 



(9) 



for reactor v c — e (this work) and accelerator u e — e [H, Q 
experiments, respectively. Accordingly, the studies of re- 
actor v e — e are expected to improve on the sensitivities of 
sin 2 #vv and (gv, 9a) at the same experimental accuracies 
as those from v e — e measurements. The relative strength 



of the three components normalized to R eX pt = 1 are in 
the ratios of 



(Rcc '■ Rnc ■ Ri 



{ (0.77 : 0.92 : 0.69) for D c - e 
' 1 (1.77 : 0.16 : 0.93) for v e - e 

_ (io) 

The stronger NC component in v c ~ e scattering is the 
physical basis of the sensitivity enhancement in the 
derivation of sin 2 6*w- 

The SM was tested and sin 2 6*w was precisely measured 
in the high energy (Q 2 > GeV 2 ) region with accelerator 
experiments on e + e~, polarized ep and vN deep inelas- 
tic processes, and in the low energy (Q 2 < 10~ 6 GeV 2 ) 
region with measurements on atomic parity violation 0] ■ 
Among them, the sin 2 6*w derived from the NuTeV ex- 
periment on vN deep inelastic scattering [TBJ was 3a 
higher than SM prediction, though the interpretations 
were complicated by strong interaction effects [2j. De- 
structive interference according to SM prediction has 
been demonstrated by accelerator v c — e scattering exper- 
iments H, Q. 

The objective of this work is to bridge the Q 2 gap in 
probing SM electroweak physics with reactor v e — e inter- 
actions. In particular, the interference effects are studied 
in this unique system. This would complement the preci- 
sion data obtained at accelerator at higher Q 2 . The mea- 
surements would place constraints on various anomalous 
neutrino interactions such as those discussed in the next 
section. 



B. Neutrino Electromagnetic Properties 

The neutrino electromagnetic interactions fl7j provide 
natural extensions to SM. The relevant parameters are 
( r2 > c ) [3 j usually called the "neutrino charge radius", 
and neutrino magnetic moments (/i„) [l8j which describe 
possible neutrino interactions with matter via the ex- 
change of virtual photons without and with the change 
of its helicity, respectively. 

Interpretations of (r 2 ) remain controversial. A 
straight-forward definition has been shown to be gauge- 
dependent and hence (r 2 J is unphysical [Jj|. However, 
there are recent attempts to define a physical observ- 
able with (r|) which give a predicted value of 
(r 2 ) = 0.4 x 10~ 32 cm 2 with the SM framework. We 
adopt in this article the more general interpretation that 
(r 2 ^ ) parametrizes contributions to non-standard interac- 
tions in neutrino scattering [2l| . 

Changes to the SM cross-sections due to (r 2 ) can be 
obtained from Eq. [5] via the replacement [l3j : 



sm 2 6*w — s- sin 0\y + ( 



V2- 



em \ / 2 



3Gf 



(11) 



where a em is the fine structure constant. 

Contributions of pi v can be described by an additional 
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Kuo-Sheng Nuclear Power Station : Reactor Building 



Reactor Pressure Vessel 




Auxiliary 
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Primary 
Containment 



Kuo-Sheng Reac )r 
Neutrino Labo: 



Suppression Pool 



FIG. 2: Schematic layout of the Kuo-Sheng Neutrino Labo- 
ratory together with the reactor core and building. 



term to Eqs. EMU 

da 
dT 



1 - T/E v 



T 



(12) 



The SM prediction of /i„ for massive Dirac neutrinos 
is extremely small (3.2 x 10~ 19 /ib where {1b is the 
Bohr magneton). However, various models with Ma- 
jorana neutrinos can give rise to [i v at the range of 
(1CP 10 — 10 -13 ) /is relevant to experiments and astro- 
physics (l8j . The most sensitive direct laboratory lim- 
its on /i„ come from high-purity germanium detectors 
at about 10 keV threshold with reactor v D [l|| [22j. At 
this low recoil energy, the \x v contributions at the present 
limit are orders of magnitude larger than those due to SM 
v c — e cross-sections. 



III. EXPERIMENTAL SET-UP 



A. Laboratory Facilities and Neutrino Flux 



displayed in Figure [3] It has been demonstrated through 
v e — proton measurements that the integrated z> -flux for 
E u > 1.8 MeV H3 and P c -spectra for E v > 3 MeV [H 
agreed with calculations to better than <3% and <5%, 
respectively. 




6 7 

Energy (MeV) 



FIG. 3: Total u c spectrum at typical power reactor operation. 



The laboratory is equipped with a 50-ton shielding 
structure depicted schematically in Figure HI consisting 
of, from outside in, 2.5 cm thick plastic scintillator panels 
with photo-multiplier tubes (PMTs) readout for cosmic- 
ray veto, 15 cm of lead, 5 cm of stainless steel support 
structures, 25 cm of boron- loaded polyethylene and 5 cm 
of Oxygen Free High Conductivity (OFHC) copper. The 
inner target volume with a dimension of 100x80x75 cm 3 
allows different detectors for various physics topics to be 
placed. Data were taken with a CsI(Tl) scintillating crys- 
tal array during data acquisition (DAQ) periods II— V. 
Each period consisted of both reactor ON and OFF data 
taking. 



A research program on low energy neutrino physics [23[ 
is being pursued by the TEXONO Collaboration at the 
Kuo-Sheng Neutrino Laboratory (KSNL), which is lo- 
cated at a distance of 28 m from Core #1 of the Kuo- 
Sheng Nuclear Power Station in Taiwan. A schematic di- 
agram is depicted in Figure The site is at the ground 
floor of the reactor building at a depth of 10 m below 
ground level, with an overburden of about 30 meter- 
water-equi valence. The nominal thermal power out- 
put is 2.9 GW. The standard operation includes about 
18 months of Reactor ON periods separated by 50 days 
of Reactor outage OFF periods when typically one-third 
of the fuel elements are replaced. 

A summary of the key information on the four data 
taking periods reported in this article is given in Table [TT1 
The evaluation of the reactor neutrino flux and spectra 
was discussed in details in Refs. [HI, [24]]. The average in- 
flux at KSNL is 6.4 x 10 12 cm _2 s _1 . A typical spectrum is 



100 cm 



75 cm 
depth = 80 cm 
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Scintillator (2.5 cm) 



Lead (15 cm) 



Stainless Steel 
Frame (5 cm) 
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OFHC Copper (5 cm) 



Movable Trolley 



V 



FIG. 4: The shielding design of KSNL. Similar structures ap- 
ply to the back and front walls. Detectors and inner shieldings 
were placed in the inner target volume. 
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TABLE II: Summary of the key information of the four data taking periods. The period numbering follows the same scheme 
as in Ref [jj] . 





Data Taking 


Reactor ON 


Reactor OFF 


DAQ 


DAQ 


Average v c 


Fiducial 


Period 


Calender Time 


Live Time 


Live Time 


Live Time 


Threshold 


Flux 


Mass 






(days) 


(days) 


(%) 


(keV) 


(10 12 cm- 2 s _1 ) 


(kg) 


II 


Feb. 2003 - Oct. 2003 


95.2 


48.4 


88.8 


100 


6.27 


43.5 


111 


Sept. 2004 - Oct. 2005 


192 


36.6 


93.4 


500 


6.50 


40.5 


IV 


Mar. 2006 - May 2007 


204.9 


43.5 


88.0 


500 


6.44 


51 


V 


June 2007 - Feb. 2008 


132.8 


27.6 


91.9 


500 


6.29 


57 


Combined 


Feb. 2003 - Feb. 2008 


624.9 


156.1 


90.4 




6.39 





Cross-Sectional View 




FIG. 5: Schematic drawing of the CsI(Tl) scintillating crystal 
array. Light output is recorded by PMTs at both ends. 



B. Conceptual Design and Motivations 

The merits of scintillating crystal detectors in generic 
low background low energy experiments have been dis- 
cussed [23]. This experiment adopted CsI(Tl) crystal 
scintillator packed in a compact array as both target and 
detector. A schematic layout is given in Figure[5l Several 
detector characteristics and design features were incorpo- 
rated [28[ which contributed to the improvement in the 
experimental sensitivities: 

(i) Proton-Free Target Region: The CsI(Tl) crystal 
is only weakly hygroscopic and does not require a 
hermetic container to seal it from ambient humidity 
(in contrast to Nal(Tl) crystal). The crystal is also 
mechanically stable and self-supporting. Therefore, 
the target region was made up almost entirely of 
CsI(Tl) (equal amount of Cs and I, with 0.15% ad- 
mixture of Tl). The other materials were the small 
amount of teflon wrapping sheets, made up of C 



and F and contributing to only about 0.13% by 
mass. There were no protons, such that possible 
neutrino-induced background from P c — p was elim- 
inated. The cross-section of this interaction is >10 2 
times higher than that of v e — c. This background 
could not be suppressed with Reactor ON/OFF 
comparisons, and could be a potential problem with 
the Savannah River experiment Q where plastic 
scintillators were adopted as target. 

(ii) Completely Active Fiducial Volume: The ab- 

sence of detector housing allowed a fiducial volume 
which was totally active. The probability of back- 
ground events to be completely measured was en- 
hanced, and this was beneficial to background un- 
derstanding and suppression. 

(iii) Complete Three-Dimensional Reconstruction: 

Each CsI(Tl) crystal module consisted of a 
hexagonal-shaped cross-section with 2 cm side 
and a length of 40 cm, giving a modular mass of 
1.87 kg. Scintillation photons were read out by 
PMTs at both ends. The sum and difference of 
the two signals provided the energy and position 
information, respectively. A three-dimensional 
reconstruction of the events was achieved. These 
information greatly enhanced the capabilities of 
background diagnostics and evaluation. In partic- 
ular, background induced by ambient radiations 
was suppressed by rejecting events at the outer 
modules or close to the PMTs. The high atomic 
number for Cs and I (Z=51 and 53, respectively) 
allowed efficient attenuation and therefore compact 
detector geometry. 

(iv) Large Mass and Expandable Detector: This 
experiment was based on a modular CsI(Tl) 
crystal array with a total mass of 187 kg. Such 
detector approach with similar target mass scale 
was also adopted in cold dark matter searches 
in the KIMS experiment [2!)j. The design can 
be easily expanded to ton-scale experiments and 
beyond. 

(v) Pulse Shape Discrimination: The light emission 

profiles of CsI(Tl) offered excellent pulse shape dis- 
crimination (PSD) between 7/e events from those 
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due to a-particles and nuclear recoils [3(1 HH . This 
allowed precise measurements of the internal con- 
taminations for background suppression and diag- 
nostics. 

(vi) Focus at High Energy Events: The reactor v P - 
spectra below 2 MeV has large uncertainties [24| . 
while ambient background dominate below the nat- 
ural radioactivity end-point of 2.6 MeV. Accord- 
ingly, only events with T > 3 MeV were studied as 
potential v e — e candidates. The low energy events 
were still recorded and analyzed for the purposes 
of calibrations and background diagnostics. 



a °- 2 
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— e/y 






1 J 




a / recoils 




i . , , i . . , i-.--.-vvv 


-vrvvvivvvrv. . 



Time (|is) 



FIG. 6: Averaged pulse shapes due to 5.4 MeV a-particles and 
7-rays of 662 keV. Nuclear recoils, as measured with scatter- 
ing with a neutron beam, give rise the same pulse shapes as 
the a-events [30l |. 



C. Detector Construction and Readout 



As depicted in Figure El the scintillating CsI(Tl) crys- 
tal detector modules were packed into a matrix array, 
with minimal inactive dead space due to the teflon wrap- 
ping sheets. The configurations varied between the dif- 
ferent DAQ periods, but the operation conditions were 
kept uniform and stable within one period. Therefore, 
each DAQ period can be taken as an independent exper- 
iment. At the end of data taking, a 12 x 9 array was 
deployed giving a total mass of 187 kg. Fiducial volume 
was defined to be the inner crystals with a separation of 
>4 cm from the PMTs at both ends. The fiducial masses 
for individual periods are given in Table [TTJ 

There were two types of crystal modules [30j from two 
production batches: (a) single crystal with 40 cm length 
were used as target placed in the central region, while (b) 
two pieces of 20 cm long crystals optically glued together 
were placed in the outer layers as active veto. The light 
output was read out at both ends of the crystal modules 
by custom designed PMTs with low-activity glass and 
diameter of 29 mm. The target array was housed inside 
a OFHC copper box of thickness 2.5 mm. Additional 
copper shielding blocks were placed on top of the box 
to fill up the inner target volume of Figure [4] The box 
was flushed with dry nitrogen to purge the radioactive 
radon gas. The CsI(Tl) array shared the target volume 
and the downstream DAQ systems with germanium de- 
tectors for magnetic moment studies and dark matter 
searches [32j . 

The electronics and DAQ systems were described in 
Ref. [33|. The DAQ system was VME-based running on 
LINUX operating system. The PMT signals were fed 
to custom-built shaping amplifiers whose output were 
recorded by Flash Analog-to-Digital Converter (FADC) 
modules at a clock rate of 20 MHz and 8-bit dynamic 
range. The DAQ trigger was generated by discrimina- 
tor set at threshold of 100 keV for P-II and 500 keV for 
P-III,IV,V, much lower than the relevant signal region. 
Signals from all sub-dominant channels with energy de- 
positions > 10 keV, as well as the PMT signals from the 
veto-panel system and various control parameters, were 
also recorded. 



A special feature of the DAQ system was the recording 
of any events delayed as much as 500 /is after the initial 
trigger. The delay time was measured with 1 [is resolu- 
tion. This allowed measurement of delayed-coincidcncc 
events due to internal radioactivity, which in turn was 
crucial to background diagnostics and suppression. The 
DAQ output was zero-suppressed, such that only those 
CsI(Tl) channels having signals within (-5 /xs, 500 /xs) 
relative to the trigger instant were recorded. 

The DAQ dead time was accurately measured by the 
random trigger (RT) events generated at 0.1 Hz uncorre- 
cted with the rest of the hardware. The typical trigger 
rate for the CsI(Tl) array was 20—30 Hz, corresponding 
to 8—12% of DAQ dead time. Data were taken with the 
germanium system in parallel in Periods II— IV, but the 
additional contributions to the DAQ dead time were only 
minor. 



D. Detector Performance 

The intrinsic performance of the CsI(Tl) crystal mod- 
ules were discussed in details in Ref. [3(J . The energy and 
position resolutions on individual module measured with 
a 137 Cs 7-source at 662 keV were 4% and <1 cm root- 
mean-square (RMS), respectively. The averaged pulse 
shapes for j/e events and a-particles in CsI(Tl) are dis- 
played in Figure [5] Separation of j/a at >99% was 
achieved by PSD down to 100 keV electron-equivalence 
energy. 

The FADC has a hardware dynamic range of only 8-bit. 
Software algorithm was devised to correct the saturated 
pulse shapes [34[ . The effective range was extended by 4 
more bits without affecting the performance parameters 
like energy resolution and PSD. The CsI(Tl) output for 
the current measurements typically saturated at about 
2 MeV, so that the events with energy <10 MeV relevant 
to the analysis were well-reconstructed. 
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A - Q R saturated only 
□ - Q L saturated only 
* - (Ql'Qr) saturated 
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FIG. 7: Typical Q L versus Q R distribution for Hl(CRV) 
events showing the background events of natural sources. Dif- 
ferent colors denote whether the PMT signals are saturated 
at their FADC readout or not. Additional software routines 
were devised to provide correct energy information for satu- 
rated events. 



IV. DATA ANALYSIS 
A. Light Collection 

The raw input to subsequent analysis were the light 
output (denoted by Ql and Qr) derived by summing 
the pedestal-subtracted FADC signals from the PMTs 
on both ends of the CsI(Tl) modules. Depicted in Fig- 
ure [7] is a typical normalized Ql versus Qr distribution. 
The selected events were those having signals only in one 
crystal (HI), with "cosmic-ray veto" (CRV) imposed and 
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software correction applied to the saturated pulses 
The different color schemes denote the status on pulse 
saturation of the two PMTs. 

Three bands along the increasing energy axis arc 
conspicuous, corresponding to background due to 7- 
rays from 137 Cs (662 keV), 40 K (1461 keV) and 208 T1 
(2614 keV). These lines were important for in situ cal- 
ibration as well as background diagnostics. The sharp 
reduction of background beyond the 208 Tl energy — the 
signal region of this measurement — is very distinct. The 
enhanced event rates at both edges indicate that most 
background sources were external to the detector. 



B. Event Reconstruction 

The objectives of event reconstruction were to provide 
measurements on energy (E) and longitudinal position 
(Z) using Ql and Qr. The calibration procedures were 



FIG. 8: The longitudinal Z-position distributions for events 
at energy corresponding to 7-lines of 137 Cs (solid histogram) 
and 40 K (dotted histogram). 



performed with in situ data, typically once every week. 

The longitudinal Z-position for the i th crystal module 
is given by: 



Z cx 



ft -Qr 



ft 



(13) 



where ft's are parameters to absorb the residual dif- 
ference in response between the left and right readout. 
The values of ft's were obtained by requiring that the 
662 keV 7-line from Cs background must be uniformly 
distributed along the length of the crystals. The pro- 
portional constants were derived by fixing the two edges 
of the distributions to be at cm and 40 cm. Typi- 
cal Z-position distributions for events at 662 keV and 
1461 keV evaluated through this prescription are shown 
in Figure The RMS resolution is 1.3 cm at 3 McV 
electron-equivalence, based on measurements with a-a 
cascade events [H| . 

The energy is described by: 



E = 



bi 



VQlx Qi 



(14) 



The parameters o^'s take into account possible differ- 
ences in the attenuation of light transmission along both 
directions, and were fixed by requiring the derived values 
of E for the 7-lines were constant and independent of Z. 
The calibration constants (a*, bi) were evaluated by a lin- 
ear fit to the 7-lines. The reconstructed energy spectra 
are depicted in Figure HI indicating RMS resolutions of 
5.8%, 5.2% and 4.0% at 137 Cs, 40 K and 208 T1 7-peaks, 
respectively. A scatter plot of the reconstructed (E,Z) 
values for a typical crystal is shown in Figure [10] The 
reconstructed energy of the various bands matched well 
to the corresponding 7-lines within the fiducial volume 
(4 cm < Z < 36 cm) . 
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Energy (keV) 

FIG. 9: The measured energy spectra at various stages of 
the analysis showing the effects of successive selection cuts. 
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FIG. 10: Energy versus Z-position scatter plot of recon- 
structed events. The prominent 7-lines matching the various 
bands are identified. 



C. Event Selection 

Neutrino-induced interactions like P e — e are of ex- 
tremely small cross-section and therefore manifest them- 
selves as "single-hit" (HI) events in only one crystal mod- 
ule uncorrelated to the rest of the system. The HI events 
were selected from raw data through selection criteria 
with CRV, "anti-Compton" multi-hit veto (MHV), and 



TABLE III: Summary of the 


suppression and si 


gnal efficiency 


factors of successive selection 


cuts within the 3— 


8 MeV energy 


range. 






Event Selection 


Background 


Signal 




Suppression 


Efficiency 


Raw Data 


1.0 


1.0 


Cosmic Ray Veto (CRV) 


0.06 


0.93 


Multi-Hit Veto (MHV) 


0.16 


0.99 


Pulse Shape 






Discrimination (PSD) 


0.34 


> 0.99 


Z-position Cut 


0.36 


0.80 



Combined 0.0011 0.77 



PSD. The CRV and MHV suppressed cosmic-induced 
background and multiple Compton scattering events 
from ambient 7-rays, respectively. The selected sample 
is denoted by HI (CRV) in subsequent discussions. The a 
and delay-cascade events from internal radioactivity [35j j , 
as well as convoluted events in accidental coincidence, 
were identified by PSD. To minimize background due to 
ambient 7-rays, an internal fiducial volume was defined. 
Events at the outermost layers of the crystal array were 
rejected, and a Z-position cut of 4 cm from both ends 
was applied to the target (inner) crystals. 

The various parameters in the calibration and selec- 
tion procedures were measured before the detector were 
assembled on site. Typically, about 10% of the in situ 
data samples uniformly distributed within a DAQ period 
were used to provide the fine adjustments. Once ob- 
tained, the optimal parameters were applied universally 
to the rest of the data set. The energy spectra at the suc- 
cessive stages of candidate event selection are depicted in 
Figure |9] 

A summary of the background suppression and signal 
efficiency factors of the cuts in the energy range of inter- 
est (3—8 MeV) are summarized in Table lIHl The signal 
efficiencies were derived from the survival fractions of RT 
events for the CRV and MHV cuts, and with the multi- 
hit Compton events for the PSD cut. The Z-position ef- 
ficiency corresponds to a 4 cm cut at both ends and were 
accounted for in the definition of the fiducial volume. 



V. BACKGROUND 

The candidate event selection procedures of Sec- 
tion ITVCl resulted in a signal-to-background ratio of 
about 1/30 at 3 MeV. The information on multiplicity, 
energy, position, cascade event timing and a/7 identifi- 
cation available for every event allowed the residual back- 
ground to be understood, analyzed and suppressed. In 
addition, the Reactor ON/OFF comparisons provided an 
independent handle to the background. These measure- 
ments were combined to improve the background evalu- 
ation which in turn enhanced the experimental sensitiv- 
ities. 
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A. Background Understanding and Diagnostics 

Several diagnostic tools which contribute to the quan- 
titative understanding of the background arc discussed in 
this section. For completeness, all prominent background 
channels are presented, though many of those are below 
the physics analysis threshold of 3 McV. 



1. Intrinsic Background 

Measurements of intrinsic radiopurity in the CsI(Tl) 
crystal with in situ data were discussed in details in 
Ref. [H. 

The isotope 137 Cs is produced artificially as fission 
waste from power reactors and atomic weapon tests. Ce- 
sium salts are soluble and can easily contaminate the raw 
materials which produce Csl. The 137 Cs contaminations 
was measured to be (1.7 ±0.3) x 10~ 17 g/g, and were 
uniform across the length of the crystals, as depicted in 
Figure [5] 

The cascade events provided measurements on 
the naturally-occurring 238 U, 232 Th and 235 U 
series, which were (0.82 ± 0.02) x 10~ 12 g/g, 
(2.23 ± 0.06) x 10~ 12 g/g and < 4.9 x 10~ 14 g/g, 
respectively, assuming secular equilibrium. The (3- 
decays of 208 T1 followed by 7's in coincidence could in 
principle lead to background in the signal region. From 
the measured level of 228 Th in the target, the contribu- 
tion of this background at 3—5 MeV was evaluated to 
be only < 11% of the expected v c — e signals. 

In addition, trace admixtures of the fission daugh- 
ter 129 I and of the naturally-occurring 40 K in the raw 
Csl powder were measured with accelerator mass spec- 
trometry techniques to be < 1.3 x 10~ 13 g/g and < 
2 x 10 -10 g/g, respectively [36] . Neutron capture on 
133 Cs produced 134 Cs at the level of - 5 x 10~ 20 g/g, as 
measured with the in situ two-hit background discussed 
in Section IVA"4l 



2. Ambient Radioactivity 



Tl Single Escape 




Multi-Hit 



Compton 



Edge 



'Cs 



40 K Energy (keV) 208 T j 



FIG. 11: Scatter plot of H2 events after cosmic-ray veto in 
0-3 MeV, showing bands on 60 Co, 134 Cs and 208 T1 single- 
escape, as well as the correlated 7's from 208 T1. 



3. Cosmic Ray Tagging Efficiency 

The cosmic-ray tagging efficiency (e M ) is the probabil- 
ity that the cosmic-ray induced events actually produce 
a "cosmic-ray tag" (CRT). The inefficiency (1 — e M ) was 
due to incomplete geometrical coverage and the light col- 
lection deficiencies of the large-area scintillator panels. 

High energy events above the end-point of (n,7) energy 
scale of about 8 MeV are all cosmic-ray induced. These 
events provided a clean sample for to be measured. 
For improved experimental control, three-hit events (H3) 
between 8—14 MeV were selected, and e M is given by 



H3(CRT) 
H3 (Total) 



(15) 



The e M averaged over all the DAQ periods was measured 
to be 92%, while the variations among periods were less 
than 1%. 



The HI spectra of Figure [9] show several 7-lines, the 
most prominent ones were those from 137 Cs, 40 K and 
208 Tl. The background dropped by several orders of 
magnitude beyond the natural radioactivity end-point 
of 2.6 McV. The cut-off at 8 MeV corresponded to the 
end-point of 7-rays emissions following neutron capture. 
The lines are crucial for energy calibration, system stabil- 
ity monitoring, and background diagnostics. Apart from 

Cs which is an intrinsic radioactivity, the other sources 
are external to the CsI(Tl) target. Distributions of the 
Z-position were heavily attenuated from the edge of the 
crystals, as illustrated for the case of 40 K in Figure [5] 



4- Diagnostics of Two-Hit Background 

Multi-hit events were unrelated to neutrino interac- 
tions and therefore provided unambiguous diagnostics to 
the background sources. Displayed in Figure ITT1 is a scat- 
ter plot of two-hit (H2) events after CRV cut. Several 
features were noted which revealed the nature and loca- 
tions of the dominant background sources, discussed as 
follows: 

(i) 208 T1 Induced Pair Production: The single es- 
cape peak following pair production of the 208 T1 
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(ii) 



2614 keV 7's can be identified. The Z-position dis- 
tribution of these events confirmed that the sources 
were external to the target. As discussed in Sec- 
tion IV Bl pair production events are crucial for 
background evaluation because of their distinctive 
topologies. 

60 Co Contaminations: It was established that the 
dominant reactor-induced radioactivity in KSNL 
was 60 Co which existed as dust in the laboratory 
area [l6| . Their contributions varied between DAQ 
periods due to different levels of contaminations 
during hardware installation. Events due to the 
correlated 7's at energy 1173 keV and 1332 keV 
from 60 Co can be located in Figure [Til They were 
uniformly distributed along the Z-position, signify- 
ing that some 60 Co dust got into the target volume 
between crystals during installation. The measured 
contamination level is 3x 10~ 20 kg -1 . However, the 
total energy of the 60 Co lines is below the 3—5 MeV 
signal region relevant to this measurement. 



(iii) Neutron Capture Induced 134 Cs: Trace 
amount of 134 Cs (r 1/2 = 2.05 yr ; Q = 2.06 MeV) 
was produced by neutron capture on 133 Cs within 
the CsI(Tl) target. It decays via /3-emission to- 
gether with two 7's of energy 605 keV and 796 keV 
in coincidence. These events were tagged in the 
H2 plot in Figure [TT1 The intensity distribution is 
uniform over the length of the crystals, verifying 
the sources were internal. The measured contami- 
nation level is 5 x 10 -20 g/g. The Q- value is below 
the physics analysis threshold and hence these 
decays would not contribute to the background of 
this measurement. 

(iv) Cascade 7-rays from 208 T1: Decays of 208 T1 are 
characterized by several 7-rays emitted in cascade. 
Coincidences of 7-rays at 510, 583 and 860 keV with 
the prominent line of 2614 keV can be identified in 
the H2 scatter plot of Figure QTJ The evaluation 
of the contributions of this channel to HI events is 
crucial to background suppression, and is addressed 
in Section EBU 

(v) Neutron Capture on 63 Cu: The main shielding 

materials in the vicinity of the target were OFHC 
copper. Neutron (n,7) capture on 63 Cu has rela- 
tively large cross-section (4.5 b), giving rise to high 
energy 7's at 7637 keV and 7916 keV. These were 
observed in HI spectrum shown in the inset of Fig- 
ure E 



5. Pair Production Event Samples 

Pair production background manifested themselves 
mostly as three- hit events (H3pp). They were selected 
by requiring two crystals each having 511 keV of energy 
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FIG. 12: Three-hit pair production spectra for CRV events, 
and further differentiating into the cosmic-ray related and un- 
related components. After scaling with the known pair pro- 
duction cross-sections, these spectra provide the in situ 7-ray 
background at the detector. 



back-to-back to the third one. These samples have dis- 
tinctive topology not contaminated by other background 
channels. Coupled with the known energy dependence of 
the pair production cross-section, these samples provided 
measurements of the in situ 7-spectra, and therefore were 
crucial for subsequent background evaluation. 

The H3pp spectrum for CRV events are displayed in 
Figure [12] This was produced by the 7-ray background 
whose contributions to the Hl(CRV) signals were eval- 
uated. There were two components to this high energy 
7-ray background: 

(i) cosmic-ray induced events with missing CR-tags, 
whose rates are given by 



H3 PP (CRV;At) = 



(1 



H3 PP (CRT) 



(16) 



where e M is the cosmic-ray tagging efficiency mea- 
sured with Eq. 1151 and 

(ii) ambient radioactivity unrelated to cosmic-rays, 
which can be evaluated with 



H3 PP (CRV; = H3 PP (CRV) - H3 PP (CRV; fi) 
also depicted in Figure [TSJ 



(17) 



B. Background Evaluation 

The experiment focused on the 3 MeV < T < 8 MeV 
energy range as the physics analysis window. The v c —e 
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FIG. 13: Measured HI spectrum and the different background 
channels evaluated in Section IV Bl The SM i> e — e contribu- 
tions are overlaid. 
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FIG. 14: Comparison of measured H2 correlated events with 
simulation results on the cascade 7-rays due to 208 Tl de- 
cays. The entry H1(CRV; Tl 7 ) denotes the predicted HI back- 
ground produced by two 7's depositing all energy only in the 
same crystal. The measured strength of the 208 T1 peak pro- 
vided consistency cross-checks to the simulations. 



signal region is expected to be at 3—5 MeV due to rapid 
decrease of the reactor P c -spcctra. 

The background diagnostics in Section IV Al demon- 
strated that convoluted 7-rays from 208 T1, cosmic-ray 
events with missing CRV tags, as well as ambient high 
energy photons could contribute to the HI background 
[Hl(BKG)]. The experimental design allowed quantita- 
tive measurement of these background which resulted in 
the extraction of the v e — c signal events with good accu- 
racy. 

The evaluation of the various background channels is 
discussed in the following sub-sections. Their contribu- 
tions are depicted in Figure [T3l where the expected SM 
v e — e spectrum is overlaid for comparison. 



1. Cascade 7 -Rays from Tl 

Decays of 208 Tl are followed by emissions of 7-rays 
in coincidence, having energy Eti(1, 2, 3, 4)= 2614.5, 
860.56, 583.2 and 510.8 keV, and at intensity ratios of 99, 
12.8, 86.2 and 25% per 208 Tl-decay, respectively. Two- 
fold coincidence manifested as H2 events were identified 
in the scatter plot of Figure [TT] Events with both 7- 
rays hitting and depositing all energy in the same crystal 
would become HI background to the v e — e signals. 

The probabilities were studied by full simulations with 
GEANT software packages [37j . incorporating realistic 
angular correlations and branching ratios for the 208 T1 
decays (38|. The sources were located at the PMTs 
and their voltage dividers, which were the only mate- 
rials other than OFHC copper and teflon in the vicinity 
of the target. 



The measured H2 distribution of E T i(l) © E T i(2, 3, 4) 
as a function of distance between the two crystals (D) 
was displayed in Figure 1141 The simulation results are 
overlaid, the normalization of which was fixed by best-fit 
to the H2 distribution for events with separation more 
than one crystal (D>1). Excellent agreement with the 
strength of the Eti(1) single-7 HI peak at 2614 keV was 
demonstrated. It served as important consistency check 
and tools for systematic studies. The data point at D=l 
denotes H2 events with hits from adjacent crystals. The 
measured intensity was significantly larger than the ex- 
pected contributions from correlated 7's due to 208 T1. 
The excess was attributed to multiple Compton scatter- 
ings at adjacent crystals from a single high energy pho- 
ton. This was reproduced in simulations studying H2 
events with single photons. 

The entry at D=0, denoted by HI (CRV; Thy), corre- 
sponds to the prediction of the HI events having the 
two 7's depositing energy exclusively in the same crys- 
tal. It was adopted for subsequent background subtrac- 
tion. The relative intensities to the 2614 keV reference 
peak, expressed as ratios of [Eti(1) © Eti(N)] /Eti(1) in 
Hl-events, are 0.13%, 0.33%, and 0.16%, for N=2,3,4, 
respectively. 



2. Cosmic-Ray Induced Background 

Once the cosmic-ray tagging efficiency (e M ) was mea- 
sured with Eg. 1151 the comic-ray induced HI background 
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with missing CR tags was derived using 



Hl(CRV ;j u) 



HI (CRT) , (18) 



similar to the H3pp channel in Eq. [T6l 



3. Ambient j-Ray Background 

This background channel H1(CRV; ]0t) is due to ambi- 
ent high energy photons emitted mostly through thermal 
neutron capture by the surrounding materials. This was 
derived directly through the H3pp (CRV; $) spectrum of 
Figure [T21 The sharp cut-off at ^8 MeV indicates the 
dominance of (n,7) processes. 

The background can be further divided into two cate- 
gories, according to the different methods of evaluation. 

(i) Full Energy Deposition (T = E 7 ): The 7-rays 
lose all its energy within a single crystal through 
multiple Compton scatterings or pair production 
with both annihilation photons fully absorbed. 
The rate is given by 



HI (CRV; ^ = 



H3 PP (CRV;^ 
H3 PP (CRT) 



HI (CRT) . (19) 



The evaluation of H3 PP (CRV;, 
Eq.HZI 



followed from 



(ii) Partial Energy Deposition (T < E 7 ): The am- 
bient 7's could undergo single Compton scattering 
after which the outgoing photons left the detector 
without further interactions. Only a fraction of the 
incident energy would be deposited in a single crys- 
tal. This background channel was studied with full- 
scale simulations using the in situ cosmic-unrelated 
H3pp(CRV; $) spectrum of Figure [T2l for flux nor- 
malization. For consistency check, the strength of 
the HI full energy (T = E 7 ) spectra of Figure [T31 
was successfully reproduced. The contribution by 
this channel to Hl(BKG) at 3—5 MeV was only 
<5% of the expected SM D c — e signals. 



4- Combined Evaluation 

It can be derived from Figure [T3] that >99% of the 
HI (CRV) events can be accounted for by the i> c ~c signals 
as well as the three dominant background channels, such 
that: 

HI (CRV) =Hl(P -c)+Hl(BKG) ; 
Hl(BKG) S (20) 
H1(CRV;T1 7 )+H1(CRV;^) + H1(CRV;^) , 

where the three contributions arc given by Fig. [T?] and 
Eqs. [T5&fT9l respectively. The sub-dominant terms in- 
clude intrinsic radiopurity and ambient 7-ray background 



TABLE IV: Summary of the sources of systematic errors 
[Say B (Source)] and their contributions to the measurement un- 
certainties [Aj, ys (£)]. The various components to the signal 
strength are summed, while those to the background subtrac- 
tion are averaged. 



Sources 


5 B ys (Source) 


A sys (0 


Signal Strength : 






$[, Evaluation 


<3% 


<0.03 


Efficiencies for Neutrino Events 


<1.3% 


<0.013 


Fiducial Target Mass 


<4% 


<0.04 


* Combined (Signal) 




<0.052 


Background Subtraction : 






Reactor OFF Measurement 


< 0.4% 


<0.06 


Background Evaluation 






0H1(CRV;T1 7 ) 


<3% 


<0.08 


©H1(CRV; ft + H1(CRV; ft 


< 1% 


<0.17 


Net 




<0.19 


* Combined (Background) 




<0.15 


Total 




<0.16 



with partial energy deposition which contributed at the 
<0.5% level of Hl(BKG). 

The H1(CRV;T1 7 ) channel was important only in the 
3—3.5 MeV energy bin. The other two channels due to 
high energy 7 interactions were dominant over the entire 
energy range of interest. Their combined contributions 
were simplified by Eqs. H6ll7ll8T S iT9l to become: 



HI (CRV; ^) + HI (CRV; 



H3 PP (CRV) 
H3 PP (CRT) 



•HI (CRT) 



(21) 

That is, the dominant contribution to Hl(BKG) was re- 
lated to the H3pp sample through a simple ratio of events 
with and without CR-tags. 



VI. SYSTEMATIC UNCERTAINTIES 

A summary of the sources of systematic errors 
[S sys (Source)] and their contributions to the measured 
£-ratio [A sys (£)] is given in Table ITVl An uncertainty of 
3% was adopted for the evaluation of the high energy 
reactor i^-spcctra. The signal efficiencies for the selec- 
tion procedures discussed in Scction FlV Cl were accurately 
measured with high statistics using the RT events. The 
fiducial mass uncertainties originated from the Z-position 
resolution of 1.3 cm. 

The systematic effects on background evaluation were 
studied with event samples unrelated to neutrino interac- 
tions accumulated over all DAQ periods. These include 
data from the Reactor OFF periods as well as those with 
energy above the 8 MeV end-point of the reactor neu- 
trino spectra. Individual methods were demonstrated to 
be able to account for the neutrino- unrelated background 
to certain accuracy levels, which were in turn assigned as 
the systematic uncertainties of those methods. 

(i) Reactor ON/OFF Comparison: The intensity of 
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TABLE V: Stability levels (8±) of the various neutrino- 
unrelated channels. 



Channels/Period 


X 2 /dof 


S± (%) 


/ua Tl Intensity : 






II 


19/13 


0.85 


III 


38/43 


0.61 


IV 


33/27 


0.81 


V 


8.4/8 


0.91 


Reactor OFF 3- 


8 MeV Hl(CRV) Rates : 




II 


15/14 


3.18 


III 


11/11 


3.51 


IV 


8.1/8 


3.60 


V 


7.2/5 


3.22 


Combined 




< 0.4% 



the 208 Tl 7- line allowed the stability of the hard- 
ware systems to be monitored and demonstrated 
to good statistical accuracies. The window within 
3—8 MeV at the Reactor OFF periods consisted ex- 
clusively of background and provided an additional 
monitor. The stability of the measured intensities 
of the 208 T1 7-line at 2614 kcV in Period III rel- 
ative to the whole-period average is illustrated in 
Figure 115b . Summary of all results are tabulated 
in Table IVl The good reduced-% 2 (x 2 /dof) indi- 
cate the data were stable within individual peri- 
ods. The hardware instability level demonstrated 
with the combined data is < 0.4%. 



(ii) 



208 T1 Induced 7-Ray Cascade Background: 

The simulation software and the normalizations 
discussed in Section IV B II were cross-checked by 
applying them to compare with the measured in- 
tensities of the 2614 keV 7-line following 208 T1 
decays for all periods, and with the HI events 
at 3-5 MeV for only Reactor OFF. The rela- 
tive deviations between the measured and pre- 
dicted rates [(R mcas - R P red)/Rmcas] are depicted 
in Figure 115b . showing consistency with zero [= 
(-0.013±0.029) at x 2 /dof = 6/11]. The fitting er- 
ror of 3% represents an upper bound of the system- 
atic uncertainties to the H1(CRV; Tl 7 ) background 
component. 

(iii) Dominant H1(CRV; /i) + H1(CRV; $) 

Background: The evaluations of the 
H1(CRV; /i) + H1(CRV; jtf) combined contributions 
in Eq. [5T] were cross-checked with measurements 
on neutrino-unrelated samples at 3—8 MeV from 
the Reactor OFF periods and at 8—12 MeV from 
both ON/OFF periods. The relative deviations 
[(Rmeas - Rpred)/R-mcas] were consistent with zero 
[= (0.0021 ± 0.0081) at x 2 /dof = 14.5/16], as 
illustrated in Figure [T5c for the combined data set. 
The fitting error of 1% represents an upper bound 
of the systematic uncertainties. 

The contributions of the individual systematic effects 
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FIG. 15: Data from which the systematic uncertainties 
were derived: (a) stability of the 208 Tl peak intensities 
in Period III; (b) comparisons of 208 T1 intensities evalu- 
ated by simulations with measured background on neutrino- 
unrelated data samples, and (c) comparisons of the evaluated 
H1(CRV; fi) + H1(CRV; //) background channels with mea- 
sured background on neutrino-unrelated data samples. 



to the 



-e cross-section measurement were then de- 



rived. The various 6 sys (Source) terms were imposed on 
the data, and the changes introduced on £ were the cor- 
responding systematic uncertainties A sys (£) listed in Ta- 
ble El 



VII. PHYSICS RESULTS 



Intermediate results of the experiment were previously 
reported [3j|. The final physics results presented in this 
section are based on the complete data set, and include 
contributions from systematic uncertainties, as well as 
improvements in the background evaluation. 
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A. Formulation 

The experimentally measured rates for neutrino events 
[Rex P t{v)\ in Eq.[7]are given by: 



RexptM = i?Hi(ON) - i?Hl(BKG) , 



(22) 



where i?ni(ON) is the measured Hl(CRV) spectra for 
Reactor ON data, and i?Hi(BKG) is the background de- 
rived from the statistical average of two different mea- 
surements: (1) Reactor OFF data, and (2) sum of the 
dominant and sub-dominant background contributions to 
Hl(BKG) in both the Reactor ON and OFF periods, the 
derivations of which are discussed in Section IVB4I 

Data from the four independent DAQ periods were 
used, combining to give a total of 29882(7369) kg-day 
of fiducial mass exposure during Reactor ON(OFF), re- 
spectively. The adopted analysis window is 3—8 MeV 
spread out uniformly over Nf, in = 10 energy bins. 



B. Cross-Section 

The cross-section ratio £ defined in Eq. [7] was derived 
with a minimum-x 2 fit, defined by 



E 



Rexptji) - £ ' RsAlji) 

A stat (i)* 



(23) 



where RsM{i) and R e xpt(j) ar e SM-expected and mea- 
sured event rates at i th bin, respectively, and A sta t(i) is 
the corresponding statistical error of the measurement. 

As cross-check, identical procedures were applied to 
the combined Reactor OFF data [i?m(OFF)], in which 
case only the predicted background was subtracted to 
provide the residual spectrum displayed in Figure 116b ,. 
Best-fit with Eg. |2"31 gave 



£(OFF) = 0.03 ± OM(stat) 



(24) 



at x 2 /dof=10.3/9, demonstrating good overall system- 
atic control of the background subtraction procedures. 

Combining all Reactor ON and OFF data from all pe- 
riods and adopting the systematic uncertainties listed in 
Table IVl the ratio 



£ = 1.08 ± 0.21(stat) ± 0.16{sys) 



(25) 



at x 2 /dof=8.7/9 was derived following Eq. [231 The mea- 
sured P e — e cross-section was consistent with the SM pre- 
diction. 

The residual and best-fit spectra are depicted in Fig- 
ure 116b . The £ ratios derived from individual peri- 
ods as well as with background subtraction by differ- 
ent methods are tabulated in Table IVII As illustrations 
using Period II Reactor ON data, the raw sample con- 
sisted of about 1.94 x 10 6 events. The analysis pro- 
cedures of Section IIV CI selected 2074 counts. A back- 
ground estimate of (2016±17±8) events was subtracted 




-0.0025 



FIG. 16: (a) The residual spectrum 

[Rhi(OFF) - R pred (OFF)] with combined Reactor 
OFF data. The best-fit to SM effects is consistent 
with £ = 0. (b) The combined residual spectrum 
[R exp t{v) = Km (ON) - .Rhi(BKG)] in the 3 - 8 MeV 
energy region. The blue and red lines correspond to the SM 
expectations and to the best-fit of the data, respectively. 



based on the various schemes in Section IV Bl result- 
ing in (57±27±8) signal events. The total i> e -e sample 
strength from all four periods combined corresponds to 
[414 ± 80(stat) ± 61(sys)\ events. 

The consistent distributions of the best-fit values and 
their errors in Table IVII together with the appropriate 
range of the x 2 /dof values indicate robustness of the anal- 
ysis procedures. These results represent a probe to SM 
at Q 2 ~ (3 x 10~ 6 ) GeV 2 and improve over those from 
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TABLE VI: Summary of the measured values of £ and xVdof 
over individual DAQ periods, as well as with the different 
background subtraction schemes in the total data set. 









€ 






X 2 /dof 


Individual Period : 














II 


1.15 


± 


0.55 


± 


0.17 


8.5/9 


III 


1.03 


± 


0.43 


± 


0.20 


8.3/9 


IV 


1.18 


± 


0.36 


± 


0.19 


7.3/9 


V 


0.97 


± 


0.42 


± 


0.20 


9.9/9 


All Periods Background Subtr 


action 












Measurement Reactor OFF 


1.25 


± 


0.43 


± 


0.08 


7.4/9 


Evaluation Reactor OFF 


1.33 


± 


0.37 


± 


0.22 


6.9/9 


Evaluation Reactor ON 


0.78 


± 


0.33 


± 


0.20 


10.3/9 


Combined 


1.08 ±0.21 ±0.16 


8.7/9 




FIG. 17: Best-fit results in (gv,gA) space and in the sin 8w 
axis from this experiment on v c — e and the LSND experiment 
on v e — e. The allowed regions are defined by their correspond- 
ing statistical uncertainties. 




FIG. 18: The measurement of interference term from best-fit 
to the data in the 3—8 MeV energy range. The solid line 
corresponds to the SM expectation of r\ = — 1. 



at x 2 /dof=8.7/9, in excellent agreement with the SM 
value of sin 2 w (SM) = 0.23867 ± 0.00016 at this low 
Q 2 (< 10 _4 GcV 2 ) range [icj. The combined uncertainty 
of ±0.039 from this measurement is less than that from 
the LSND accelerator v c — e experiment of ±0.051. The 
improvement is due to the enhancement factors favoring 
D c — c processes, as indicated in Eq. [9] 

The interference term was probed using Eq. [3] The 
best-fit value of the sign-parameter rj is 



'/ 



-0.92 ± 0.30(stat) ± 0.24(sys) 



(27) 



at x 2 /dof=8.8/9. The residual spectrum showing 
{Rexpt — Rcc — -Rnc) is displayed in Figure [HI with the 
expected spectra for 77 = 0, ±1 overlaid. The results veri- 
fied destructive interference in the SM D e — e interactions. 



previous reactor neutrino experiments 

I2H1H- 

C. Electroweak Parameters 

The constraints on the coupling constants {gvi9A) 
were derived by a minimum-^ 2 two-parameter fit on 
Eq. [3] with a formulation similar to that of Eq. [23] The 
allowed regions are presented in Figure [TT] Results from 
the accelerator experiment with v e Q are overlaid. The 
complementarity of the v c — e and D c — e processes is read- 
ily seen. 

The weak mixing angle was derived with best-fit on 
Eq. EJ giving 

sin 2 6» w = 0.251 ± 0.031(stai) ± 0.024(s?/s) (26) 



D. Neutrino Electromagnetic Properties 

To include possible effects due to \l v and (r 2 o ), the 
expression of Eq. [21] was modified to 

fl(/i„;<4» = ifei(ON)-[JZBM(i/) + JZ H1 (BKG)] • (28) 

The non-standard interaction parameter (r 2 ) as de- 
fined in Eq. [TTJ was measured to be 

(r 2 c ) = [0.61±1.30(stoi)±1.01(sys)] x 10~ 32 cm 2 (29) 

at x 2 /dof=8.7/9. This can be translated to bounds of 

- 2.1 x 10~ 32 cm 2 < (r 2 J < 3.3 x 10~ 32 cm 2 (30) 

at 90% confidence level, an improvement over the current 
limits by the LSND experiment with v e — e @: —2.97 x 
10" 32 cm 2 < (r 2 J < 4.14 x 10~ 32 cm 2 . 
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TABLE VII: Projected statistical sensitivities on £ and 
sin 2 #w under various realistically achievable improvement to 
the experiment. 

Improvement A stat (£) A stat [sin 2 fl w ] 



This Work 


0.21 


0.031 


Improved Feature : 






A. x 10 Data Strength 


0.07 


0.010 


B. Background Reduction 






Bl: >99% Cosmic-Ray Efficiency 


0.12 


0.018 


B2: x jq Reduction in 






Ambient & 208 T1 7's 


0.16 


0.024 


* With Both B1+B2 


0.05 


0.007 


All Features A+B1+B2 Combined 


0.015 


0.0022 



Similarly, the best-fit value for [ij 1 is 

M„ 2 = [0.42 ± 1.79(stat) ± lA9(sys)] x 10~ 20 fj.% (31) 

at x 2 /dof=8.7/9, which corresponds to a limit of the D e 
neutrino magnetic moment of 

fx 9e < 2.2 x 10~ 10 /i B (32) 

at 90% confidence level. This is, however, less stringent 
than the best published limit of Mp c < 0-74 x 10~ 10 Mb 
with germanium detector at 12 keV analysis thresh- 
old jig. 

VIII. SUMMARY AND PROSPECTS 

We report in this article an improved measurement 
of reactor v c with the atomic electrons at the Q 2 ~ 
10~ 6 GeV 2 range. Complementary and comparable 
sensitivities on the SM electroweak parameters were 
achieved as those measurements with accelerator v c at 
higher Q 2 . 

The detector concept allowed complete three- 
dimensional event reconstruction, with which wc 
demonstrated that the background above 3 McV could 
be identified, studied and accounted for to the level of 
~1% accuracy. The background understanding and sub- 
traction enhanced the experimental sensitivities beyond 
the conventional Reactor ON— OFF comparisons. 

The sensitivities can be further enhanced. As illus- 
trations, the projected improvement under various re- 
alistically achievable assumptions are summarized in Ta- 
ble lVHl Electromagnetic calorimeters using CsI(Tl) with 



tens of tons of mass have been constructed, such that the 
target mass is easily expandable. As shown in Figure IT3l 
the dominant background above 3 MeV were all exter- 
nal to the target scintillator. Accordingly, they will be 
attenuated effectively through self-shielding in a target 
with bigger mass. The incorporated features listed in Ta- 
ble |VIl] correspond to 10 times increase in data strength 
and >10 times suppression in background. The statisti- 
cal accuracies can be improved to 1.5% and 0.9% for £ 
and sin 2 #w, respectively. 

As indicated in Table IIV| the systematic errors on 
background subtraction are related to the actual back- 
ground level, such that they will also contribute to 
A S3/S (£) at the < 2% level under the assumption of Ta- 
ble IVIII Modest improvement on the evaluation of re- 
actor neutrino spectra will attain similar accuracy. To 
eliminate the errors in fiducial mass, active light guides 
with a different scintillating crystal can be coupled to 
both ends of the CsI(Tl) target crystal. A good candi- 
date is the pure Csl crystal. The vast difference in the 
scintillation decay times (~10 ns versus ~1000 ns for Csl 
and CsI(Tl), respectively) [3(| makes the definition of an 
inner target volume simple and exact using PSD tech- 
niques. 

The projected sensitivities of such experiments are sim- 
ilar to those estimated with a large liquid scintillator tar- 
get [4ll ]. and can complement the sin 2 #w measurements 
with the high energy accelerator experiments. In partic- 
ular, these experiments can probe the anomalous NuTeV 
results [l5| at comparable sensitivities but with a differ- 
ent neutrino interaction channel and at a low Q 2 [42T) . 
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